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ABSTRACT 

We perform a joint analysis of X-ray and Sunyaev-Zel’dovich effect data using an analytic model that describes the 
gas properties of galaxy clusters. The joint analysis allows the measurement of the cluster gas mass fraction profile 
and Hubble constant independent of cosmological parameters. Weak cosmological priors are used to calculate the 
overdensity radius within which the gas mass fractions are reported. Such an analysis can provide direct constraints 
on the evolution of the cluster gas mass fraction with redshift. We validate the model and the joint analysis on 
high signal-to-noise data from the Chandra X-ray Observatory and the Sunyaev-Zel’dovich Array for two clusters, 
A2631 and A2204. 
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1. INTRODUCTION 

Galaxy clusters trace the growth of structure in the uni- 
verse. Their abundance and evolution is critically sensitive 
to underlying cosmological parameters such as Oj/, rr x , and 
the dark energy equation of state parameter w. Recent work 
has focused on using galaxy clusters to constrain cosmol- 
ogy, including dark energy constraints from X-ray and joint 
X-ray/Sunyaev-Zel’dovich (SZ) measurements of the gas mass 
fraction (Allen et al. 2008; Ettori et al. 2009; LaRoque et al. 
2006), cosmological parameter constraints from the growth 
of structure via X-ray (Mantz et al. 2008, 2010; Vikhlinin 
et al. 2009), and SZ cluster surveys (Vanderlinde et al. 2010; 
Marriage et al. 2011; Sehgal et al. 2011; Muchovej et al. 2011; 
Williamson et al. 2011; Benson et al. 201 1). 

In this paper, we present a method for the joint analysis 
of X-ray and SZ cluster observations using a self-consistent 
analytic model for the physical properties of the intracluster 
medium (ICM; Bulbul et al. 2010). The model provides ana- 
lytic expressions for the radial density, temperature, and pres- 
sure profiles, and is therefore simultaneously applicable to both 
X-ray and SZ observables. The joint analysis allows measure- 
ment of the cluster gas mass fraction without the need to impose 
external priors on cosmological parameters such as the Hubble 
expansion rate H(z). Such an analysis applied to a sample of 
clusters can directly probe the evolution of cluster gas mass 
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fractions with redshift. We demonstrate the method using high 
signal-to-noise data from Chandra and Sunyaev-Zel’dovich Ar- 
ray (SZA) observations of two clusters, A2631 and A2204. 

The method developed in this paper combines all available 
Chandra X-ray data (both imaging and spectroscopic) with SZ 
observations, using these to determine the angular diameter dis- 
tance and cluster mass with minimal cosmological assumptions. 
We chose the Bulbul et al. (2010) model for this analysis since 
it describes the three thermodynamical cluster properties — 
density, temperature, and pressure — with a consistent set of 
parameters that are both readily interpreted. The temperature 
profile linking density and pressure is both easily calculated and 
observationally motivated. This work differs from Mroczkowski 
et al. (2009) who used SZ and X-ray imaging data only with a 
simplified, core-cut form of the Vikhlinin et al. (2006) model 
to describe the X-ray density and the Nagai et al. (2007) pa- 
rameterization of the SZ pressure, and for which the inferred 
temperature profile did not reduce to a compact, accessible 
expression. 

This paper is structured as follows: Section 2 describes the 
data reduction and analysis, Section 3 the modeling of X-ray 
and SZ data, and Section 4 the joint analysis of the X-ray and 
SZ data. We present and discuss our conclusions in Section 5. 

2. DATA REDUCTION 

We selected a non-cool core cluster (A2631) and a cool-core 
cluster (A2204) with high-quality X-ray and SZ observations to 
demonstrate the method of analysis. 
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Table 1 

Cluster Observations 


Hasler et al. 


Cluster 

z 

nu 

R.A. 

Deck 

Chandra Observations 


SZA Observations 




(10 20 cm' 2 ) 

(J2000) 

(J2000) 

ObsID 

Time 

(ks) 

On-source Time 
(hr) 

FWHM 

(arcsec) 

P.A. 

(deg) 

rms Noise 3 
(mJy) 

A2631 

0.27 

3.55 

23:37:40.1 

+00:16:33 

3248/11728 

25.0 

16.1 

152 x 117 

17.2 

0.4 

A2204 

0.15 

5.67 

16:32:47.2 

+05:34:32 

7940 

72.9 

19.6 

157 x 115 

-7.7 

0.4 


Notes. a FWHM (full width at half-maximum of the synthesized beam), P.A. (position angle of the synthesized beam), and rms noise are for short baselines 
«2U). 


Table 2 

Radio Sources in Cluster Fields 


Cluster 

Pointing Center 


30 GHz Source 


FWHM 

(arcsec) 

P.A. 

(deg) 

rms Noise b 
(mJy) 

R.A. 

(J2000) 

Deck 

(J2000) 

Source 

Aa a 

(arcsec) 

A <5 a 

(arcsec) 

Flux 

(mJy) 

A2631 

23:37:38.8 

+00:16:06.5 

i 

21.3 

36.5 

3.7 

26.5 x 16.6 

42.7 

0.25 




2 

205.0 

-130.0 

0.5 




A2204 

16:32:46.88 

+05:34:32.4 

1 

0.4 

1.2 

7.0 

21.1 x 18.4 

-82.1 

0.22 




2 

-417.8 

-360.1 

21.6 







3 

195.0 

-130.1 

0.7 





Notes. 

a A a and AS are the offsets from the pointing center. 
b FWFIM, P.A., and rms noise are for long baselines (> 2k/.). 


2.1. Chandra Imaging and Spectroscopy 

The Chandra X-ray data are in the form of event files, which 
we use to generate both images and spectra. Additional blank- 
sky composite event files are used for background subtraction. 
The event files are reduced using CIAO 4.3.1 and CALDB 4.3, 
following the reduction procedure described in Bulbul et al. 
(2010). Details of each cluster observation can be found in 
Table 1. 

The X-ray images in the 0.7-7 keV band are used to measure 
the X-ray surface brightness profile of the cluster. To subtract the 
background from the surface brightness, we rescale the blank- 
sky image to match the cluster surface brightness in a peripheral 
region that is free of cluster signal. The peripheral regions we 
chose are at a distance ^500 arcsec from the cluster center for 
A263 1 (corresponding to 2.7 Mpc in the standard flat ACDM 
cosmology) and ^550 arcsec (1.7 Mpc) for A2204. 

Spectra are extracted in annular regions centered at the peak of 
the X-ray emission. These regions cover an area out to the radius 
where the surface brightness profile reaches the background, 
which is near 1-500 (the radius within which the average density 
is 500 times the critical density at the cluster redshift) for these 
observations. From the blank- sky data, background spectra are 
also extracted and processed. We then rescale the blank-sky 
spectra to match the count rate of the cluster spectra in the 
9.5-12 keV band. In this band, Chandra has no effective area for 
the detection of photons, and the detected counts originate from a 
particle background that is time variable. Flickox & Markevitch 
(2006) showed that while the flux within the 2-7 keV and 
9.5-12 keV energy bands can vary with time, the ratio of the two 
bands remains constant. Subtracting the blank-sky data rescaled 
by the higher-energy band therefore accurately removes the 
background from the lower-energy band. 

After rescaling the blank-sky spectra and removing the back- 
ground from our cluster data, residuals may still be present in 
the soft 0.7-2 keV energy band. These soft X-ray residuals may 
be due to Galactic and extragalactic emission, and may vary as 


function of position (e.g., Snowden et al. 1997) and time (e.g., 
Takei et al. 2008). For each cluster observation, we use a pe- 
ripheral region that is free of cluster emission — the same region 
used to rescale the background images — to determine whether 
soft residuals are present after the blank-sky background has 
been subtracted. We detect the presence of soft residuals in both 
clusters. The residual spectra are fit using a phenomenological 
model that includes a power law and a plasma emission model, 
and this model is rescaled by area and included in the spectral 
fit for each annulus (e.g., Snowden et al. 1998; Nevalainen et al. 
2005; Maughan et al. 2008). 

2.2. Interferometric Observations with the SZA 

The two clusters were observed with the SZA, an eight- 
element interferometer designed for the detection and imaging 
of the SZ effect (Figure 1). Each antenna in the array is 
3.5 m in diameter and has a primary beam FWFIM of 1017 
at the center frequency of the observing band (31 GHz). For 
these observations, six antennas were closely packed together 
to provide sensitivity to arcminute-scale SZ signals, and the 
remaining two antennas were placed further out to constrain the 
flux contributions from unresolved radio sources, as described 
in Muchovej et al. (2007). The unflagged on-source time for 
A2631 was 16.1 hr and 19.6 hr for A2204. The details of the 
observations are given in Table 1, and the radio sources detected 
in each field are listed in Table 2. In the analysis of the cluster 
SZ effect described below, the parameters of the SZ decrement 
and the radio sources are fit simultaneously. 

The SZA data are reduced using a set of routines written 
in MATLAB 14 that constitute a complete pipeline for flagging, 
calibrating, and reducing visibility data. The reduction pipeline, 
described in Muchovej et al. (2007), converts the data to 
physical units and corrects for instrumental phase and amplitude 
variations. Data are flagged for corruption due to bad weather, 


14 http://www.mathworks.com/products/matlab 
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Right Ascension (J2000) Right Ascension (J2000) 


Figure 1 . SZ contours overlaid on X-ray false color images of A2631 (left) and A2204 (right). The Chandra X-ray surface brightness data are from the energy range 
0. 7-7.0 keV. The color bars reflect the number of counts detected by Chandra in the 0.7-7 keV band, with a pixel size of 1.97 arcsec. The SZ data are from the SZA, 
and the contour levels are (+2,— 2,— 4,— 6,— 8, ...) times the rms noise (see Table 1). The FWHM of the synthesized beams for these SZ observations are shown in the 
lower left comer of each image. 

(A color version of this figure is available in the online journal.) 


sources of radio interference, and other instrumental effects 
that could impact their quality. The pipeline outputs calibrated 
unflagged visibilities, i.e., components of the Fourier transform 
of the sky brightness multiplied by the primary beam response, 
along with their corresponding statistical weights and positions 
in the Fourier (u-v) plane. 

3. MODELING THE X-RAY AND SZ DATA 

The X-ray observable is spatially resolved spectroscopy, 
and temperature and metallicity of the intracluster plasma are 
measured using the X-ray spectroscopic data. The X-ray surface 
brightness is defined as 

Sx = . ,, [ n 2 e A ee (T e , A)di, (1) 

4tt(1+z) 3 J 

where i is the line of sight through the cluster, n e is the electron 
density, T e is the electron temperature, A is the metallicity, 
and A ee (T e , A) is the X-ray cooling function (in units of 
counts cm 3 s -1 ) as a function of electron temperature and 
metallicity. The density, temperature, and metallicity can vary 
along the line of sight. 

The observable from the SZ data is the amplitude of the 
spectral distortion of the cosmic microwave background (CMB) 
in the direction of the cluster. This distortion is due to inverse 
Compton scattering of CMB photons off electrons in the ICM, 
and results in a decrement in the CMB brightness temperature 
at frequencies of < 218 GHz. The magnitude of the decrement 
is proportional to the electron pressure integrated along the line 
of sight (Sunyaev & Zel’dovich 1972): 

A7sz = Tcmb [ crx f(x, T e )n e — e —dl, (2) 

J m e c A 

where 7cmb is the temperature of the CMB, fix. T e ) contains 
the frequency dependence of the SZ temperature signature using 
the relativistic corrections provided by Itoh et al. (1998) and 
Nozawa et al. (2006), cs\ is the Thomson cross section, m e is the 
electron mass, and c is the speed of light. 


Assuming spherical symmetry, the line-of-sight integration 
element dl relates to the angular element (in radians) as 
dl = D\d0, where Da is the angular diameter distance. From 
Equations (1) and (2), we find 


S x a j 

f D A n 2 e A ee (T e , A)d0 

(3) 

AT$z oc j 

f D A n e T e d6. 

(4) 


The combination of X-ray imaging spectroscopy and SZ obser- 
vations can be used to simultaneously measure the distribution 
of the electron density, the electron temperature, and the angu- 
lar diameter distance (e.g., Hughes & Birkinshaw 1998; Grego 
et al. 2000; Reese et al. 2002; Grainge et al. 2002; Saunders et al. 
2003; Bonamente et al. 2006). For further discussion on the SZ 
effect and its use for cosmology, see reviews by Birkinshaw 
(1999) and Carlstrom et al. (2002). 

We describe the density and temperature profiles of the hot 
plasma in galaxy clusters using the model proposed by Bulbul 
et al. (2010): 

n e (r) = n e0 ^(.r,r s ,P) n r~ l ol (5) 


T(r) = T 0 (p(r, r s , P)r COO h 

where 


<P(r, r s , P) 


1 (1 + r/r s )P- 2 - 1 

(P - 2) r/r s ( 1 +r/r s y~ 2 ' 


( 6 ) 

(7) 


_ a + (r/ r cool y 

^cool 1 . x .. 5 

1 + (r/r coo i) y 


( 8 ) 


n e o is the normalization of the pressure profile, n is the poly tropic 
index, r s is the scale radius, Tq is the normalization factor for the 
scaling of the temperature profile, y is the slope of the cooling 
function, and a is the cooling parameter which ranges from 0 
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Radius(arcsec) Radius(arcsec) 

Figure 2. X-ray surface brightness and temperature profiles for A2631 (left column) and A2204 (right column). Top panels: surface brightness profiles where the 
black points are derived from the X-ray images, the red line shows the best-fit model, and the green line shows the background level determined from the blank-sky 
observations. The residuals show the fractional difference between the model and the data. A 1% systematic uncertainty has been added in quadrature to the datapoints. 
We plotted the surface brightness profiles beyond the fitted region to show the agreement between the cluster emission and background. Bottom panels: temperature 
profiles where the red line shows the best-fit model and the dark, medium, and light gray regions show the 68%, 95.4%, and 99.7% confidence levels obtained from 
the model fits. A 10% systematic uncertainty has been added in quadrature to the temperature bins. 

(A color version of this figure is available in the online journal.) 


to 1 . One attractive feature of these models is that they provide 
a simple analytic form for the electron pressure: 


Pe(r)= P e0 <P( r ,r,,pT +] , (9) 

where P e o = n e okTo is the pressure normalization. The free 
parameters of the model, which are used to jointly fit the X-ray 
and SZ data, are n e o, To, r s , r coo i, a, /5, y, n, and the distance/^. 

Parameter estimation is done using a Monte Carlo Markov 
chain (MCMC) method described in Bonamente et al. (2004). 
Correlation among model parameters is a feature of most 
analytic models, including the one used in this paper. Parameter 
correlations can result in low acceptance rates and thus slow 
convergence of the Markov chains (e.g., Gilks et al. 1996). 
In our implementation, steps of the MCMC are proposed 
along the set of directions which diagonalize the covariance 
of the posterior distribution, determined via singular-value 
decomposition (SVD) of an initial test chain, resulting in 
efficient exploration of the parameter space (see the Appendix). 

3.1. X-Ray Data Analysis 

The annular bins in the temperature profile (Figure 2) were 
chosen by starting with an initial 10" bin and then increasing 
each bin by 50% of the width of the previous bin to give 
roughly the same counts per bin. Following the analysis of 
the systematic uncertainties for the X-ray data described in 


Bulbul et al. (2010), we adopt a 1% systematic uncertainty 
on the count rate of each bin of the surface brightness profile 
and a 10% systematic uncertainty on the temperature of each 
spectral region as discussed in Section 4.5.1. Figure 2 shows 
the temperature and surface brightness profiles, along with the 
best-fit models, for the fit to the X-ray data only of A2631 
and A2204. Since our model has the same parameters for both 
density and temperature, the surface brightness profile carries 
a larger weight in the fit. The model fits are acceptable for 
both clusters to within the plotted errors, which include the 
systematic uncertainties associated with the surface brightness 
and temperature discussed in Section 4.5. For A2631, there is 
insufficient signal to constrain n and f J > simultaneously. In this 
case, we fix j3 — 2, which is equivalent to assuming that the 
total mass follows a Navarro et al. (1997) profile at large radii 
(see Section 4.3). 


3.2. SZ Data Analysis 


After the removal of compact radio sources in the cluster field, 
the visibilities V v {u, v) measured by the SZA can be related 
to the Fourier domain equivalent of the integrated Compton-y 
parameter Y(u, v) (see, e.g., Mroczkowski et al. 2009). This is 
defined as 


Y{u, v ) = 


V v (u, v) 
g(x)lo ’ 


GO) 
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u— v radius(kX) u— v radius(kX) 

Figure 3. SZ visibility profiles for A2631 (top) and A2204 (bottom) plotted as a function of u—v radius (\! u 1 + v 2 ). The plots show the real components of the 
measured Y{u,v) along with the best-fit model. 


where g(x) corrects for the frequency dependence of the SZ 
flux, and 7o = 2(k bPcmb) 3 / (/*c) 2 is the primary CMB intensity 
normalization. Figure 3 shows Y(u, v) along with the best-fit 
model, for the fit to the SZ data only of the two clusters. 

4. JOINT ANALYSIS OF X-RAY AND SZ DATA 

We first perform a consistency check of the determination 
of the pressure profiles from the X-ray and SZ data. We then 
focus on determinations of the angular diameter distance and 
the radial profile of the gas mass fraction using the consistent 
parameterization of density, temperature, and pressure provided 
by the ICM model for the joint analysis of the X-ray and SZ 
observables. 

4.1. Consistency of X-Ray and SZ Measurements 
of the Electron Pressure Profiles 

X-ray and SZ observations provide independent measure- 
ments of the radial distribution of the electron pressure. The 
X-ray observables are the electron temperature and surface 
brightness, the latter depends on the square of the electron den- 
sity according to Equation (1); the SZ observable, on the other 
hand, is directly proportional to the electron pressure integrated 
along the line of sight (Equation (2)). 

The two observables can be affected by different sources 
of systematic uncertainties. For example, the presence of non- 
thermal X-ray emission (e.g., Million & Allen 2009; Bonamente 
et al. 2005; Sarazin & Lieu 1998) could result in the increase 
of the X-ray emission above the level of the thermal gas, 
and radio emission from cluster halos (e.g., Brunetti et al. 
2007) may partially fill the SZ decrement. Another source of 
systematic uncertainty is the assumption of spherical symmetry 
in the analysis (e.g., Sulkanen 1999), which would result in 
a different measurement of the pressure from X-ray and SZ 
observations. A discussion of sources of systematic uncertainty 
in the analysis of X-ray and SZ observations is presented in 
Section 4.5. A comparison of the pressure profiles from SZ and 
X-ray observations is therefore useful to determine the presence 
of sources of emission that can cause differences between the 
two measurements. 

We perform a joint fit to the X-ray and SZ data using 
Equation (9), with the normalization of the SZ pressure model 
independent of the X-ray density and temperature normaliza- 
tions. The common parameters in the models ( r s , fi>, and n ) are 


Table 3 

Pressure Normalization Values and Integrated FCrsoo) Values 2 


Cluster 

Peo(SZ) P eo (X) 

(1CV 10 erg cm -3 ) 

Peo(SZ)/Peo(X) 

A2631 

A2204 

L00±°;!! 

9.90±°:“ 

i 21+ 015 

1 - Z1=t 0.14 

9 71± 0 - 47 
' i=c 0.47 

0.82±°;[» 

l-02±gS 


Lph.szLsoo) 

(to- 

Lph,x('‘50o) 

- H ) 

Lph.szLsoo)/ T sp h,x(r50o) 

A2631 

A2204 

9-i3±l;g 

44.97±|*j 

n-i3±};|5 

43.93± 3 ; 03 

0.82±°;i 3 

i 02+ 010 
i .uziCq 09 


Note. a Statistical and Chandra calibration systematics are included in the 
measurement of masses. 


linked between the two data sets, thus requiring the X-ray and 
SZ pressure profiles to have the same shape. In this analysis, we 
adopt the angular diameter distance appropriate for the cluster 
redshift in a ACDM model with h = 0.73, Q-m = 0.27, and 
£2 A = 0.73. 

The pressure inferred from the X-ray and SZ measurements 
are within 20% of one another for both clusters, consistent with 
the statistical and systematic effects (Tables 3 and 6); results 
for a larger sample of clusters can be found in Bonamente 
et al. (2012). Note that the measurement of the ratio of X-ray 
pressure to the SZ pressure depends on the choice of the Hubble 
parameter, since the pressure normalizations are degenerate with 
the value of Da assumed in the analysis. 

4.2. Direct Measurement of the Angular Diameter Distance 

We also perform a joint X-ray and SZ analysis which enables 
us to place direct constraints on the angular diameter distance 
without using priors on the cosmological parameters (e.g., 
Hughes & Birkinshaw 1998; Birkinshaw 1999; Reese et al. 
2002; Bonamente et al. 2006). For this analysis, we link the 
shape parameters and the pressure normalizations between the 
X-ray and SZ data, and allow D A to vary. For A263 1 and A2204 
we measure D A = 798.9±^Mpc and D A = 575.3±^Mpc, 
and both values are consistent with those calculated using a 
standard ACDM cosmology at the ler level. The measurement 
of the D a for A2204 is also in agreement with that of Bonamente 
et al. (2006). The measurement of the angular diameter distance 
for a given cluster is affected by a number of systematic effects 
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Figure 4. Top panels (left A2631, right A2204): gas mass profiles determined from the joint analysis of Chandra / SZA observations for A2631 and A2204. Middle 
panels: total mass profiles. Bottom panels: gas mass fraction profiles. The dashed lines are the 68% confidence level at each radius. Gray areas show the measurements 
at radii 7*2500 and 7 * 500 , obtained by marginalization over the cosmological parameters. Chandra calibration systematics are included in the measurements. 


(Bonamente et al. 2006), and the agreement of Da with the 
ACDM value is expected for a large sample but not necessarily 
for individual clusters, as it is for the two clusters in this paper. 

4.3. Radial Profiles of the Gas Mass Fraction Independent 
of Cosmological Parameters 

By using direct measurement of Da as described in 
Section 4.2, we can also obtain radial profiles of the gas mass, 
total mass, and the gas mass fraction without the need to use 
priors on the cosmological parameters. The gas mass M gas 
is computed by integrating the gas density profile within the 
volume, 


is computed assuming hydrostatic equilibrium between the 
gravitational mass and the thermal pressure of the gas: 


M tot (r) = 


^Ttp-,r] 

03-2) 



1/(1 -fl-r/rA 
(1 +r/r s )P~ l ) 


tcool(f), 

( 12 ) 


where p,- = (Tok(n + 1)(/J — l))/(4jtGpm p rj). 

Figure 4 shows the radial profiles of the gas mass, total mass, 
and gas mass fraction for A2631 and A2204. The uncertainties 
reflect the fact that Da is also measured directly from the data, 
and that no assumption about the value of the cosmological 
parameters 7 /q, or was made. 


Mgas = 4tt fi e m p f n e (r)r 2 dr — 4jt jl e m p Da f n e (6)6 2 d9, 

( 11 ) 

where p is the mean molecular weight (calculated assuming 
metal abundances of 0.3 solar; Anders & Grevesse 1989), m p 
is the proton mass, and dr = D A d6. The total mass M tot 


4.4. Measurement of the Gas Mass Fraction 
at an Overdensity Radius 

In cosmological applications (e.g., via the distribution of / gas 
with redshift; Allen et al. 2008) / gas is typically measured within 
an overdensity radius. The radius r A is defined as the radius 
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Table 4 

Results from Joint X-Ray /SZ Analysis 


Cluster 




Model Parameters 






n e o 

r s 

n 


V 

r cool 

a 

y b 

Da 


(KT 2 cm -3 ) 

(arcsec) 



(keV) 

(arcsec) 



(Mpc) 

A2631 

0.78±»;! 2 

26 1 1+74.6 
zoi . 13=50 2 

9.6± 2 j 

2.0 

7.6±{; 2 




799±267 

A2204 

3.90±°; 23 

22.7±}J 

6-9±l:l 

i 27+0-10 
1 /=I= 0.08 

14.9±JJ 

20.0±qj 

0.17 ±0.01 

2.0 

575± 47 





Cluster Masses 0 







Masses Evaluated at A = 2500 



Masses Evaluated at A = 500 



Cluster 


A^gas 

Mtot 

/gas 

r A 

Afgas 

M tot 

/gas 



n 

(10 13 m o ) 

(to 14 M 0 ) 


(") 

(10 13 Mq) 

(10 14 m o ) 



A2631 

103.2±!^ 

1 71 + 240 

i 37 _l1.16 

1 •* 5 /=c 0.65 

0 124+ 0081 
U.lZH-3=oo60 

289.3± 44 ;° 

8.35±“g 7 

6.16±j;M 

0 131± 0103 
u. id i3=o 065 


A2204 

23 1 .3±g3 

5 03± 108 

D.U33=i 0? 

4 09± 0 39 

H-.uy3= 034 

0 122± 0019 
U. 1ZZ3=0 018 

492.6± 23 ; 3 

13-09 ± 2 2 f 4 

7 95± 1-20 
'■ yj ~ n 0 .92 

0 16 i+o° 26 

u. 1013=0 021 



Notes. 

a The reader is cautioned that To is not a global temperature, but rather a model parameter in Equation (6). 
b The parameter y is fixed in the model. 

c Statistical and Chandra calibration systematics are included in the measurements. 


within which the average matter density of the cluster is A times 
the critical density of the universe at the cluster’s redshift: 



M tot (r a) 

yAftb)' 


( 13 ) 


where p c (z) = {3H^E 2 {z))/{%nG) is the critical density of the 
universe, Hq is the Hubble constant, and E 2 {z)— 12 vi ( 1 + z) 3 + 
Q A + Qk(l + z) 2 in the ACDM model. 

The joint X-ray and SZ analysis provides cosmology- 
independent constraints on Da and on the radial profile of /g as (r) 
(Sections 4.2 and 4.3). The radius r A and therefore all quantities 
calculated out to this radius retain a cosmological dependence 
through the factor p,(z.) appearing in Equation (13). In the fol- 
lowing, we describe a method to marginalize the measurement 
°f /gas(i"A) over the cosmological parameters, which results in 
a weak cosmology dependence of our joint measurements of 

,/gas (^Al- 
in the following, we adopt the standard flat Friedman- 
Robertson-Walker (FRW) cosmological model, parameterized 
by H a and 12 V i , with 12 A = 1 — 12 M , which is known to provide 
a good fit to other cosmic distance data (Allen et al. 2008; 
Freedman et al. 2009; Percival et al. 2010). Within this model, 
a given D A value corresponds to a curve in the (Hq, Q m ) plane 
described by 


and the low redshifts in question, this procedure corresponds, 
to very good approximation, to the use of uniform priors on 
Hq and 12 m- We note, however, that this correspondence does 
not necessarily hold for higher redshifts (z ~ 1) or for other 
cosmological models. The measurements of masses and / gas at 
7*2500 and r $ oo are indicated as gray areas in Figure 4 and are 
listed in Table 4. 

For comparison, we also provide the measurements of masses 
and / gas using the X-ray data with a fixed D A calculated from 
the standard flat ACDM cosmology, and priors on 12 m from 
the Wilkinson Microwave Anisotropy Probe measurements of 
Komatsu et al. (2011). The results are reported in Table 5. 

4.5. Sources of Systematic Uncertainty 

The sources of systematic uncertainty on the gas mass frac- 
tion and the integrated X-ray and SZ pressure are listed in Ta- 
ble 6. The individual errors are added in quadrature to determine 
the total systematic uncertainty on / gas and the integrated pres- 
sure. Systematics from Chandra instrument calibration and the 
X-ray background are included in the fitting of the data, and the 
value of masses and / gas in Tables 4 and 5 account for these sys- 
tematics. Since this joint analysis method can be applied to larger 
samples of clusters, we also indicate whether the impact of each 
source of systematic uncertainty is reduced with sample size. 

4.5.1. Instrument Calibration 


D A (l+z)J 0 [ 12 m ((1 + f ) 3 — 1 ) — 1] 1/2 ’ 

where spatial flatness is assumed. For each step in the Markov 
chain of each cluster, we use the corresponding value of D A 
to generate a consistent pair of cosmological parameters by 
drawing a random value of Qm from a uniform distribution 
on [0,1], and calculating the associated value of Hq from 
Equation (14). We then use these parameters when calculating 
r A and the other fit parameters of the chain. This approach 
has the advantage of using minimal prior information on the 
cosmological parameters, with £2 m free to range between 0 
and 1, and D A measured directly from the data. This method 
properly accounts for the covariance between r A , / gas (r A ), D A , 
and the other parameters of the fit. For the flat FRW model 


Chandra’s ACIS effective area has a spatially dependent 
non-uniformity at the level of ±1% and therefore we add a 
±1% uncertainty to the surface brightness data. We also adopt a 
±10% uncertainty on the temperature measurements to account 
for uncertainty in the low-energy calibration of the effective area 
(see, e.g., Bulbul et al. 2010). These uncertainty estimates for 
the Chandra calibration are folded into the mass measurements 
reported in Tables 3-5. 

Frequent observations of Mars are used to calibrate the SZA 
absolute flux scale; we employ the Rudy (1987) flux model 
which has an estimated absolute calibration uncertainty of ±5%. 
The stability of the instrumental gain is ±3%, as determined 
from repeated calibrator measurements in SZA survey fields 
(Muchovej et al. 2011). The absolute calibration and instrumen- 
tal gain yield a global ±6% uncertainty on the SZA calibration. 
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Table 5 

Results from X-Ray Analysis 


Hasler et al. 






Model Parameters 





Cluster 

n e 0 

r s 

n 

p 

7b a 

r cool 

a 

r b 

£>A b 


( 1 (T 2 cm -3 ) 

(arcsec) 



(keV) 

(arcsec) 



(Mpc) 

A2631 

0 7r+ 003 

u - /r=c 0.03 

248.0±«™ 

9.27±};f 3 

2.0 

8.3 ±«; 2 




840.0 

A2204 

4 12 ± 0-23 
^• az=1= 0.25 

22.5 ±}; 2 

6.76±‘;» 

1 ao_i_0.06 
1.3o ± 0 06 

1 5 . 0 ± q ! 

20 . 0±° : 2 

f) 1'7_|_0.08 
UA /=t 0.09 

2.0 

526.0 





Cluster Masses 0 







Masses Evaluated at A = 2500 



Masses Evaluated at A = 500 



Cluster 

7*A 

4^gas 

4^tot 

/gas 

7*A 

^fgas 

Mtot 

/gas 



(") 

(10 13 M e ) 

(10 14 Mo) 


(") 

(10 13 M g ) 

(10 14 m o ) 



A2631 

113.9±®;® 

2 45± 0 - 26 
z.hou=027 

1.98±!J; 32 

0 124± 0010 
U - 1Z ^ =C 0.008 

309.9± 2 ° : 2 

10.34±° : 3 ° 

8 00± L33 

0 129=t 0019 
uMzynioois 


A2204 

234.4± 3 J 

4 12 ± 009 
^• az=c 0.09 

3.75±«;» 

0 1 1 o ± 0003 

u. i 1 o=n 0 003 

499.0±|‘j 

10.64±«;g 

7 27± 0-53 
'• z /= l = 0.53 

0 146± 0007 
u. ih-o=h 0 007 



Notes. 

a The reader is cautioned that To is not a global temperature, but rather a model parameter in Equation ( 6 ). 
b Parameters y and Da are fixed in the model. 

c Statistical and Chandra calibration systematics are included in the measurement of masses. 


Table 6 

Sources of Uncertainty 


Source 


7*2500 



>*500 



Effect on / gas 

Effect on Pressure 

Effect on / gas 

Effect on Pressure 


(%) 

SZ (%) 

X-ray (%) 

(%) 

SZ (%) 

X-ray (%) 

Kinetic SZ effect 

±4 

±4 


±4 

±4 


Radio point sources 

±2 

±1 


±2 

±1 


Asphericity 

±20 

±10 

±10 

±20 

±10 

±10 

X-ray background 

±5 


±1 

±9 


±2 

SZA calibration 

±10 

±6 


±10 

±6 


Hydrostatic equilibrium 3 

-9 



-11 



Model assumptions 

±10 

±3 

±3 

±10 

±5 

±5 

Helium sedimentation 

+ 10 


-4 

+5 


-2 

Total systematic 

±28 

±13 

±11 

.27 

“ C 29 

±14 

±12 

Chandra calibration uncertainties 15 
Surface brightness 

±10 






Temperature 

±1 







Notes. 

a Uncertainty is theoretically motivated by Lau et al. (2009) as discussed in Section 4.3.4. 

b These systematic uncertainties are added to the data prior to the fit, and their effect is included in the derived masses and pressure at 
all radii. 


We rescaled the SZA data by the ±6% uncertainty on the SZA 
calibration and compared the measurements to the original anal- 
ysis. We found a ±6% systematic uncertainty on the pressure 
and ±10% systematic uncertainty on the gas mass fraction at 
7*2500 and 7 * 500 . The uncertainty associated with instrument cali- 
bration does not average down with sample size. 

4.5.2. Kinetic SZ Effect 

Reese et al. (2002) report that for a cluster with a temperature 
of 8.0 keV and with a typical velocity along the line of sight of 
300 km s -1 (Watkins 1997; Colberg et al. 2000) the kinetic SZ 
effect would be ±4% of the thermal SZ for 30 GHz observations. 
Accordingly we use a ±4% uncertainty due to the kinetic SZ 
effect on the gas mass fraction and the SZ pressure profiles 
for measurements at 7*2500 and 7*500- This source of uncertainty 
averages down by a factor of the square root of the size of the 
sample. 


4.5.3. Radio Source Contamination 

Undetected radio sources not accounted in the modeling could 
lead to a biased measurement of the SZ decrement. We use the 
Faint Images of the Radio Sky at Twenty-centimeters (FIRST) 
database as a reference for locating compact radio sources within 
10' of the cluster center. Most radio sources that will affect the 
30 GHz data (rms noise ~ 0.25 mJy) will have counterparts in 
the FIRST survey (rms noise of 0.15 mJy at 1.4 GHz); inverted 
spectrum sources may not have counterparts at 1 .4 GHz and will 
affect our measurement of the SZ decrement, but they comprise 
a small fraction of the source population (Muchovej et al. 2010, 
Figure 3). 

We determine the effect of undetected point sources by 
placing a radio source model at each FIRST source, fixing 
the position and marginalizing over the flux. In the pressure 
model (Equation (7)), we fixed the parameters r s , n, ff and Da, 
and let P eo be free. We compare the pressure profiles with the 
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original analysis (see Table 2) and find a ^ 1% difference in the 
pressure of each cluster. Therefore we conservatively apply a 
1% uncertainty on the pressure and 2% uncertainty on the gas 
mass fraction at nsoo and rsoo. 

4.5.4. Asphericity 

Although we assume a spherical model in our analysis, most 
clusters do not appear to be circular in shape in X-ray or 
radio observations. LaRoque et al. (2006) report a 10%-20% 
uncertainty in the measurement of the gas mass fraction due to 
asphericity; therefore we use ±20% as a conservative estimate 
for measurements at rysoo and rsoo ■ This uncertainty also 
averages down by a factor of the square root of the sample 
size, as shown in Sulkanen (1999), provided the selection of the 
sample is unbiased with respect to cluster shape. 

4.5.5. Hydrostatic Equilibrium Assumption 

The assumption of hydrostatic equilibrium at large radii 
results in an underestimate of the total mass. This is due to the 
presence of non-thermal pressure which can bias hydrostatic 
equilibrium measurements of the total mass. According to Lau 
et al. (2009) the total mass of a relaxed cluster such as A2204 
will be biased by —6% at r 25 oo and —8% at rsoo, and for 
unrelaxed systems such as A263 1 by —9% at r 2 soo and — 1 1 % at 
rsoo- Therefore, we adopt a systematic uncertainty of —9% and 
— 11% in our error analysis for measurements at r 2 soo and rsoo, 
respectively. 

4.5.6. X-Ray Background 

The X-ray background is determined from the AC1S blank- 
sky composite event file. We normalize the blank-sky back- 
ground level for each observation using an emission-free region 
on the ACIS detector. We adjusted the background normaliza- 
tion factor by a factor of ±2<r and propagated this through the 
analysis, and found that this produces a ±2% uncertainty on 
the background count rate. This uncertainty affects the surface 
brightness and temperature measurements resulting in a ±2% 
and ±9% uncertainty on the gas mass fraction measurements at 
r 25 oo and rsoo, respectively, and a ±2% uncertainty on the X-ray 
pressure profiles at both radii. This uncertainty averages down 
by a factor of the square root of the sample size. 

4.5.7. Systematics Associated with the Use of 
the Bulbul et al. (2010) Models 

The Bulbul et al. (2010) model assumes a polytropic relation- 
ship between the ICM density and temperature at large radii. To 
estimate uncertainties associated with the polytropic assump- 
tion, we compare our X-ray masses for the two clusters in our 
sample with those calculated using the Vikhlinin et al. (2006) 
model, which provides an independent parameterization of the 
thermodynamic quantities. From this comparison, we find that 
the gas mass fraction measurements varies by ± 10% at all radii 
between the Bulbul et al. (2010) and Vikhlinin et al. (2006) mod- 
els. We estimate the uncertainty on the corresponding pressure 
profiles by comparing the integrated pressures between the two 
models, and the comparison results in an uncertainty of ±3% 
at r 25 oo and ±5% at rsoo- We consider these uncertainties as a 
rough estimate of the systematics associated with the Bulbul 
et al. (2010) model. 


4.5.8. Helium Sedimentation 

The effect of helium sedimentation may be an additional 
source of systematic uncertainty. In our measurements we 
assume that the hydrogen to helium ratio is uniform throughout 
the cluster. However, theoretical studies (Fabian & Pringle 1977; 
Rephaeli 1978) suggest helium sedimentation effects may affect 
cluster mass measurements. Peng & Nagai (2009) find that 
the bias in gas mass fraction from the presence of helium 
sedimentation is less than 10% at ~r 2 soo and negligible at rsoo. 
Accordingly, we estimate a systematic uncertainty of 10% at 
T 2500 an d ^5% at rsoo for the measurement of the gas mass 
fraction. Bulbul et al. (2011) applied the Peng & Nagai (2009) 
helium sedimentation simulation model to a sample of clusters 
and demonstrated the effects on the gas mass and total mass. The 
integrated pressure is proportional to the gas mass, and we use 
the values from Bulbul et al. (2011) to determine upper limits 
to the systematic uncertainty in the measurement of pressure of 
-4% at r 250 o and -2% at r 500 . 

5. CONCLUSIONS 

We demonstrate the use of the Bulbul et al. (2010) cluster 
model for simultaneous fitting of X-ray data and SZ effect data. 
The model employs a compact parameterization that relates the 
three primary thermodynamic quantities by the ideal gas law at 
all radii. We consider X-ray data from Chandra and 30 GHz 
SZ data from the SZA for both clusters, A263 1 and A2204, and 
find that the model adequately captures the radial variation in 
both the X-ray surface brightness and SZ Compton-y profiles. 
For all clusters, separate determinations of the electron pressure 
from the X-ray and SZ data yield profiles that are statistically 
consistent. 

Joint analysis of the X-ray and SZ data provides a direct 
measure of D A , the angular diameter distance to the cluster, that 
is independent of cosmology. For both clusters, this analysis 
yields a measure of D A that is consistent with the standard 
ACDM values at the \a level. Using the measured angular 
diameter distance as a constraint between Hq and Qm, we 
marginalize over the implicit cosmology dependence of the 
overdensity radius to obtain estimates of / gas at risoo and rsoo 
that are only weakly dependent on £2 m- 

We discuss possible sources of systematic errors in the / gas 
determination and find that most will be mitigated if / gas is 
averaged over a large sample of clusters. A sample spanning 
a large redshift range can be used to constrain the evolution 
of / gas with redshift, and for constraining cosmological models 
with clusters (e.g., Sasaki 1996; Pen 1997; Reese et al. 2002; 
Allen et al. 2004, 2008; LaRoque et al. 2006; Bonamente et al. 
2006; Ettori et al. 2009). 
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the NSF under a cooperative agreement, and by the CARMA 
partner universities. S.M. acknowledges support from an NSF 
Astronomy and Astrophysics Fellowship; C.G. and S.M. from 
NSF Graduate Research Fellowships; and D.P.M. from the 
NASA Hubble Fellowship grant HF-5 1259.01. Support for this 
work was provided for T.M. by NASA through the Einstein 
Fellowship Program, grant PF0-1 10077. 
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n 

Figure 5. Scatter plot for parameters of the Bulbul et al. (2010) model applied to the X-ray analysis of A2631. Note the strong correlation among some of the 
parameters, especially n and r s . For clarity, only every 100th parameter in the chain is plotted. 


APPENDIX 

MCMC REPARAMETERIZATION USING A 
SINGULAR- VALUE DECOMPOSITION METHOD 

Correlation among model parameters is a common fea- 
ture of analytic models such as the beta model (Cavaliere & 
Fusco-Femiano 1976), the Vikhlinin et al. (2006) model, and 
the Bulbul et al. (2010) model used in this paper. The MCMC 
method for the analysis of X-ray and SZ effect data described 
in Bonamente et al. (2004) accounts for this correlation, and 
therefore correlation is not an issue when evaluating integrated 
quantities such as masses and Y values, and their uncertainties. 
Strong parameter correlation, however, may cause the MCMC 
to be inefficient in its sampling of parameter space (see, e.g., 
Gilks et al. 1996, page 90), requiring long chains with low ac- 
ceptance rate because of the poor mixing. A common solution 
is the use of an SVD (e.g.. Press et al. 1992) to perform a lin- 
ear transformation of the parameters to reduce the correlation 
among model parameters, and increase of the rate of accep- 
tance in the MCMC. For the X-ray analysis of the Chandra 
data of A2631 shown in Table 5, the four model parameters 
(n e o, r v , n, and fi) are transformed into four SVD parame- 
ters (svdo through .vvc/ 3 ), and the usual Metropolis-Hastings 
MCMC is applied to the SVD parameters. The accepted param- 


Table 7 

Correlation Coefficients for the X-Ray Analysis of A263 1 


Parameter 

n eo 

r s 

n 

Txo 

^eo 


-0.85 

-0.82 

0.31 

r s 



1.00 

-0.37 

n 




-0.35 

T x0 





Parameter 

svdo 

svdi 

svd\ 

svd-$ 

svdo 


-0.01 

-0.14 

0.21 

svd\ 



0.02 

-0.04 

svd2 




-0.15 

svd3 






eters are then transformed back to the original Bulbul et al. 
( 2010 ) model parameters for which we calculate integrated 
quantities. 

The effect of the reparameterization is shown in Figures 5 
and 6 . The strong correlation present between certain pairs 
of parameters, especially r s and n, is absent from the SVD 
parameters. 

The values of the correlation coefficients for the original 
parameters and the SVD parameters are shown in Table 7. 
With this reparameterization, we obtain an acceptance rate of 
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Figure 6. Scatter plot for the SVD parameters of the same chain as in Figure 5. 


SVD2 


approximately 30%, which is a factor of few higher than the 
typical acceptance rate obtained using the original parameters. 
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